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ABSTRACT 

Two types of nova evolutions can be realized in low-mass white dwarfs of ^ 0.5-0.7 Mq, i.e., an evolution 
with optically thick winds like in usual classical novae, or an another type of evolution without them like in 
the symbiotic nova PU Vul. The latter type is characterized by spectra of no indication of strong winds as 
well as a long-lasted flat optical peak in its light curve. We propose a transition from no-optically thick wind 
evolution to usual evolution with optically thick winds as a new outburst model for slow novae that show a 
relatively long-lasted multipeak phase followed by a wind phase like in the slow novae V723 Cas, HR Del, and 
V5558 Sgr We calculated nova envelopes with one-dimensional approximation of the companion's effects and 
found that when the companion star is deeply embedded in the extended nova envelope, the structure of static 
envelope approaches that of the optically thick wind solution. Thus, the transition from static to wind solution 
is triggered by the effect of the companion. The transition occurs in a close binary nova like V723 Cas, but is 
not triggered in a long period binary like PU Vul. We reconfirm our previous results that the frictional energy 
deposition is negligibly small in almost all of hydrogen/helium novae because of the low envelope density at 
the orbit. 

Subject headings: binaries: close — novae, cataclysmic variables — stars: individual {W723 Cas, HR Del, 
V5558 Sgr) — stars: winds, outflows 



1. INTRODUCTION 

Nova is a thermonuclear runaway event on a mass-accreting 
white dwarf (WD). During a nova outburst, the envelope on 
the WD expands considerably and engulfs its mass-donating 
companion unless it is a wide binary. As the companion orbits 
in the nova envelope, frictional processes between the enve- 
lope and the companion produces thermal energy. In 1980's, it 
was believed that these frictional processes play an important 
role to eject the envelope and t o shorten the nova duration by 
a factor of 10 or mor e (MacD onald et al. II1985I : lOvio et al. I 
II990t IShanker et al. Iil99 1). This idea was originally pro- 
posed to reconcile observed short durations of classical novae 
1 yr) with long timescales of nuclear burning 1000 yr). 
However, the drag energy deposition is not effective for mass 
ejection in the common envelope phase of novae because the 
envelope masses are too s mall to produce drag -energy enough 
to eject the envelope (Kat o & Hachisul[l991allbh . 

In the beginning of 1990's, opacity tables were revised 
(OPAL opacity: Iglesias & Rogers 1996), of which a promi- 
nent peak at log T{K)^ 5.2 is strong enough to accelerate op- 
tically thick winds even in relatively less massive WDs of 
~ 0.6 Mq. This wind is as massive as 10""* Mq yr"' or 
more. The nova duration is drastically reduced so that the 
theoretical duration be comes com parable to the observed one 
(e.g.. lKato & Hachisuil994lPriainik & Kovetz 1 11995). Once 
the wind occurs, frictional effects are ineffective, because the 
density at the orbit is too low to produce large energy de- 
position. Moreover, the wind is accelerated deep inside the 
orbit and the wind velocity has akeady reached the escape 
velocity at the orbit. Thus, once the optically thick winds oc- 
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cur, the presence of a companion has virtually no effects on 
the envelope structure nor on the duration of nova outbursts 
( lKato&HachisullT99 4). 

There are observational indications that frictional process 
is not effective in nova outbursts. If the frictional process is 
very effective in mass ejection, there should be a clear depen- 
dence of nova light curve shapes on the orbital period, but no 
such properties are known. For example, the recurrent nova 
RS Oph shows a ver y fast development of optical light curve 
jHachisuetal. 1 12006) similar to that of U Sco. RS Oph is 
a long period binary of PqA = 456 days jFekel et al. 2000), 
so that the WD photosphere does n ot reach the co mpan- 
ion even at the maximum expansion jHachisuetal. I I2OO61) . 
Therefore, the frictional processes of companion do not work 
in the whole period of the nova outburst. On the other 
hand, U Sco is a short period binary of Porb = 1-23 days 
(Schaefer& Ringwald 1995) and the companion star is en- 
gulfed deeply in the envelope during the outburst. The 1999 
and 2010 outbursts of U Sco were densely observed, but there 
is no indication for additional acceleration nor enhanced mass 
ejection due to the companion motion (Hachisu et al. 2000). 

We have analyzed a number of nova light curves with 
the universal decli ne law based on the optically thick winds 
(Hachisu & Kat ol |20M 12 007. 2010; Hachisu et al. 200^ 
f200&: ,Kato & Hachisu lll994. 2009: .Kato et aL.2009.) . If the 
companion motion strongly contributes to the acceleration of 
matter, there should be a change in the decline rates of light 
curves when the companion reappears from the extended en- 
velope of the WD. However, none of them show a particu- 
lar change in the light curves before and after the epoch that 
the companion reappeared from the envelope. We refer two 
examples of classical novae, with a densely observed light 
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curve and known orbital period: V15 00 Cyg (Pprb = 0-14 days , 
see the light curves in Figure 13 of iHachisu & Katoll2006h : 
the companion emerged on days ~ 50 but nothing had hap- 
pe ned, and also V458 V ul (forb = 0.59 days, in Figure 26 
of IHachisu &Tatol l20lO^: the companion emerged on days 
~ 60 but nothing happened. These are good counterevidence 
against the effectiveness of frictional process. Therefore, we 
emphasize that frictional processes are not effective in novae 
whenever strong optically thick winds blow. 

In our theoretical models, optically thick winds always oc- 
cur in nova outbursts on massive WDs (> 0.7 Mq) but do not 
occur in less massive WDs (< 0.5 Mq). Between them, i.e., 
~0.5-0.7 Mq, both types of nova evolutions, with/without 
optically thick winds, are realized for the same WD mass 
dKato & Hachisu Il2009l) . When no optically thick winds are 
accelerated, nova outbursts show a slow evolution with a long- 
lasted flat optical maximum, in contrast with usual classical 
novae that show a sharp optical maximum cau sed by a rapid 
evolution due to wind mass loss (lKatoll20"Tl1) . The first ex- 
ample of such no-optically thick wind evolution is the sym- 
biotic nova PU Vul, in which a flat optical peak lasted 3000 
days and showed very quiet spectra indicating no strong winds 
dKato et al. 2 0ll|) 

Kat o & Hachisu I (|2069) presented an idea that a nova out- 
burst started in a state of no-optically thick winds (we call this 
"static evolution") and then possibly changes to be in a state 
of optically thick winds ("wind evolution"). They also sug- 
gested that such a transition may accompany violent activities 
such as oscillatory behaviors in the light curves of some slow 
novae. 

Here, we present an idea that a companion star plays a 
role in changing the nova envelope structure from "static" to 
"wind." If this occurs in low-mass WDs, it may relate to the 
peculiar light curves of slow novae. As an application of such 
a transition to slow novae like V723 Cas, we have examined 
effects of a companion star with one-dimensional (ID) ap- 
proximation of nova envelopes. In Section|2] we briefly intro- 
duce our method. Numerical results are presented in Section 
[3] Applications to slow novae appear in Section|4] Discussion 
and summary follow in Sections |5]and|6] 

2. ENVELOPE MODEL 

We consider low-mass WDs in which no optically thick 
winds are accelerated in the beginning of nova outbursts. This 
may happen in low-mass WDs (< 0.7 Mq) under some con- 
ditions of the initial envelope mass and chemical composi- 
tion (see Kato & Hachisu 2009, for more details). After a 
thermonuclear runaway sets in, the envelope greatly expands 
and the optical brightness increases to reach a flat maximum. 
We approximate the nova evolution by a sequence of hydro- 
static solutions when no optically thick winds occur, and by a 
sequence of steady-state solutions with optically thick winds 
when it occurs. In these sequences of solutions the envelope 
mass is decreasing due to nuclear burning and wind mass loss 
if it occurs, from which we calculate evolution timescale. This 
quasi-evolution is a good approximation because the ther- 
mal/dynamical timescale between two successive solutions 
are short enough compared with the evolution timescale of 
novae. 

For static solutions we solve the equations of hydrostatic 
balance, mass continuity, radiative diffusion, and conserva- 
tion of energy, from the bottom of the hydrogen-rich enve- 
lope through the photosphere assuming spherical symmetry. 
For wind solutions, we solve the equation of motion assuming 



steady -state, instead of hydrostatic balance (i Kato & Hachisu I 
119941) . Convective energy transport is calculated in static so- 
lutions using the mixing length theory with t he ratio of mixing 
length to pressure scale hight a = 1 .5 (see IKato & Hachisu I 
2009, for effects of the a parameter on static solutions). The 
occurrence of opti cally thick w inds is detected by the con- 
dition described in iKato I (Il985l) . In our numerical calcula- 
tions, we adopt more than 8000 meshes and always insert 
~ 1 700 meshes in the interacting region between r = Ro^b - Rn 
and r = Rorh + Ra, where /?orb is the position of the companion 
from the WD center and R.^ is the accretion radius as defined 
below. The chemical composition of the envelope is simply 
assumed to be uniform with the composition of X = 0.55 , Y = 
0.23, CNO = 0.2, and Z = 0.02, which are representative val- 
ues for slow novae (see Table 1 of H achisu & Kato 120061) . We 
assume an 0.55 Mq WD and an 0.4 Mq main-sequence (MS) 
companion that fills its Roche lobe unless otherwise stated. 
This 0.4 Mq is close to an upper limit of the companion mass, 
0.44 Mq, in a binary system with thermally stable mass trans- 
fer {q = Mcomp/A^WD 0-8)- We assume such a large mass 
because a more massive companion has larger effects on the 
WD envelope. 

Three effects of the companion on the WD envelope are 
considered here: (1) spun-up by the companion motion, (2) 
gravity of the companion star, and (3) drag luminosity due to 
frictional energy deposition. We incorporated these effects in 
our computer code as explained in the following subsections. 

2.1. Centrifugal Force 

In the very beginning of a nova outburst, the envelope has 
so small specific angular momentum compared with the spe- 
cific orbital angular momentum. After the envelope expands 
to a giant size, it begins to counterrotate on the rotating frame 
with the companion due to conservation of the local angu- 
lar momentum. The outer part of the envelope may be spun 
up due to frictional effects of the companion motion. We in- 
corporate the effect of centrifugal force in our ID code. We 
modified the equation of hydrostatic balance as 



1 dP 
p dr 



GMr , 
^ + ruj\ 



(1) 



where cu is the angular velocity. This is a ID approximation 
that may represent envelope structure in the equatorial plane. 
Here we introduce the effective mass parameter / and rewrite 
Equation ([TJ as 

I'J^-^f. (2) 



p dr 



GM, 



(3) 



The factor / represents the degree how much the gravity is ef- 
fectively reduced due to the centrifugal force. In other words, 
the effective mass decreases to fM^. We, hereafter, call / the 
mass reducing factor 

We assumed a simplified rotation law for the envelope; a 
rigid rotation in the inner part of the envelope (r < Rn) and 
a Keplerian-like rotation in the outer part (r > Rq). Since 
the envelope is spun up by the companion motion inside the 
outer Lagrange point, we define the rigid rotation radius as 
Rn = Rorb+ 1.2 X Rcomp bccausc the outer Lagrange point is 
located at ^ 1 .2 x /?comp from the companion center. 

For the inner part (r < Rq), we assume a solid rotation with 
an angular velocity of ojq = ijVoTb/Roib, here 77 is a parameter 
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that represents the ratio of the angular velocity of the enve- 
lope to that of the companion, fl = Vorb/^orb- In other words, 
7] means the ratio of the envelope velocity to the companion 
velocity at the orbit. 

For the outer part (r > Rq) we assume that the rotation ve- 
locity of the envelope varies as r~'/^ as the Keplerian velocity 
does. In this case, the angular velocity uj changes as r"^/^, and 
can be expressed as a; = rjVorbRorbRli^'''^^^- Using Kepler's 



law V^u ■ 

orb 

can be rewritten as 



: G(MwD+A^comp)/^orb, the Specific centrifugal force 



Then / becomes 



"orb 



Rn G(MwD+Meon,p) 



f=l-v'i^r, for r>Rn. 

"orb 



(4) 



(5) 



Here we use M,. =MwD+-^comp at r > Rq. Thus, / is constant 
at r > Rn- Approximations of these centrifugal forc e and ID 
approximation are discussed in Sections [5 .21 and |53] 

2.2. Gravity of the Secondary 

In our ID approximation, we have simply assumed that 
the companion is distributed uniformly within a shell of r = 
Roih i^comp- Since the envelope mass (10-^ to 10"*' M©) IS 
much smaller than that of the companion (Mcomp = 0.4 Mq), 
the change of M,- is essentially due to the mass distribu- 
tion of the companion, i.e., Mr = Mwd inside the radius of 

= ^orb -^comp, and Mr = MwD +A^comp outsidc the radius of 

= ^orb +^comD- These tre atments are essentially the same as 
those in lKato&Hachisu 'l 1994). 

Figure [1] shows examples of the angular velocity u, mass 
distribution Mr, mass reducing factor /, and effective mass 
/Ml for our envelope models. Here, / is unity in the deep 
inside the envelope (r < Rorb), but 0.971, 0.741, 0.578 and 
0.418 in the outer region (r > Rn) for r/ = 0.1, 0.3, 0.383, 
and 0.45, respectively. When 77 = 0.383 (/ = 0.578), we have 
(0.55 Mq + 0.4 Mq) X 0.578 = 0.55 Mq at r> Rn, i.e., the 
companion gravity is almost canceled by the centrifugal force. 
A small peak of the effective mass appears near logr (cm)^ 
11.0 because Mr is constant (i.e., Mwo + A^comp) outside r = 
^orb + ^comp, whcrcas u! begins to decrease from a bit outside, 
at Rn = Roib + 1 .2 x Rcomp 

2.3. Drag Luminosity 

We treat the drag luminosity in spherical (ID) ap- 
proximation in__die__saine__way as in the other previ- 
ous works (iTaam & Bodenhei mer I [T989t iLivio et al. 1 119901: 
iKato & Hachisu 11199 1a b, 1991). The drag luminosity gen- 
erated in the region from r to r+Sr can be approximated as 



^idrag = /5(Vorb-Ve) SS, 



(6) 



where p is the density of the envelope, Vorb - Ve is the relative 
velocity between the envelope and the companion star, 5S is 
the cross-sectional area between a circular strip from rto r+Sr 
and a circle with the radius 7?a, the center of which is located 
at the orbit, r = /?orb- Here the modified accretion radius, R-^, 
is defined as 



R.= 



Ro 



l + (/?o/2//)2' 



(7) 
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Fig. 1 . — Adopted rotation law \ogu] (green solid lines), gravity reducing 
factor / (black dotted lines), mass distribution Mr (red solid line), and effec- 
tive mass, fMr (blue solid lines), for a binary consisting of an 0.55 Mq WD 
and an 0.4 Mq MS companion with r;=0.1, 0.3, 0.383, and 0.45, from upper 
to lower for / and fMr, but from lower to upper for \oguj. The correspond- 
ing envelope structures are shown in Figure|4] The left/right edge of each line 
corresponds to the bottom/photosphere of the envelope. The location of the 
companion (logr (cm)=l 0.869) and its size are indicated by a small open cir- 
cle with a short horizontal bai'. The orange solid and dash-three-dotted lines 
denote a tria l rot ation law and the corresponding effective mass, respectively 
(see Section l53t . 

where // is th e local density scale height, and TJq is the gener- 
alized |BondD (Il952h radius defined by 



Ro = 



IGMr., 



[(Vorb-Ve)2 + C2]' 



(8) 



here Cj is the sound speed. In our ID calculation the drag 
luminosity, (5Ldrag, is re-distributed over the entire spherical 
shell between r and r+ 6r. The total drag luminosity is 



^drac 



P(V< 



orb 



e)' 5S, 



(9) 



where the integral region is from Rg^ - /?a to /?orb + ^a ■ We take 
the companion radius itself, Rmmp, as the accretion radius in- 
stead of /?a when R^ is smaller than the companion radius. If 
both the density and /?a are constant in space, the drag lumi- 
nosity becomes the standard expression, i.e.. 



Ldrag = T^RlpiVorb-Ve)^ 

dTaam & Bodenheimer|[T989l:lLivio et al. Ill990h . 



(10) 



3. NOVA ENVELOPES WITH A COMPANION STAR 

3.1. Structure of Static Envelopes 

Figure|2]shows distribution of the temperature, density, dif- 
fusive luminosity, and the local Eddington luminosity. We 
do not include the effects of a companion for blue and green 
lines. The local Eddington luminosity (dotted lines) is defined 

L^-*-^, (1.) 
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Fig. 2. — Distributions of the temperature T, density p, diffusive l umi - 
nosity Lr, and the Eddington luminosity L^dd defined by Equation lilt . 
Each line corresponds to the envelope with the photospheric temperature of 
logrph(K) = 3.9 (thick black lines), 4.0 (middle thin black lines), 4.3 (left 
thin black Hues), and 4.6 (red lines). The location of the companion and 
its size are indicated by the small open circle with a short horizontal bar 
For comparison, four solutions without the companion effects are added: 
log rph(K) = 3.81 (green thick lines: right) and 4.3, 4.6 (green thin lines: right 
and left, respectively), and 4.9 (blue thin lines). The two envelopes of thick 
hues (with/without the companion effects) have the same envelope mass of 
4.1 X 10-^ Mq. 

where k is the o pacity. We used the OPAL opacity 
(llglesias & Rogers1[l9 96). As the opacity is a function of the 
temperature and density, and Mr is also the function of the ra- 
dius (see Figure [U, the Eddington luminosity is also a local 
variable. 

The blue and three green lines in Figure |2] depict the en- 
velope solutions of logTph (K)=4.9, 4.6, 4.3 and 3.81 for 
the case of no companion star. These envelopes have essen- 
tially the s ame structure exc e pt for the surface region. As 
reported by iKato & Hachisu I (l2009h . a static envelope has a 
core-halo structure of the density in which a large density 
inversion layer develops deep inside the photosphere. The 
density-inversion region corresponds to the super-Eddington 
region. The local Eddington luminosity has the deepest local 
minimum between logr (cm) ~ 10.0 and 11.0, correspond- 
ing to the peak of the OPAL opacity at logT (K) ^ 5.2. This 
density-inversion arises in order to keep hydrostatic balance 
in the super-Eddington region {Lr > iEdd) as expected from 
the equation of hydrostatic balance (Kato & Hachisu 200^. 
Inefficient convections occur in the region of Lr > Ledd but are 
unable to carry all of the diffusive energy flux. 

Black and red lines in Figure|2]indicate the models in which 
the companion's gravity and drag luminosity are incorporated. 
The three envelopes denoted by the black lines have very simi- 
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Fig . 3 . — Same as Figure[2]but for static solutions with different companion 
masses. The photospheric temperature is logTph (K)= 3.85 for all the solu- 
tions. Mcomp= 0.0 (blue), 0.1 (brown), 0.2 (pink), 0.3 (green), and 0.4 Mq 
(black sohd). Diffusive luminosity for Mcomp= 0.2 and 0.3 Mq models are 
omitted to simplify the figure. The location of the companion are indicated 
by the small open circles with a short horizontal bar for Mcorap= 0. 1, 0.2, 0.3, 
and 0.4 Mq from left to right. Dashed and dotted lines indicate static solu- 
tions with enhanced drag luminosity by a factor of 2.8 and 10, respectively 
(see Section lS^ for details). 

lar structures at log r (cm) < 11.5. The drag energy deposition 
is small as explained below. Therefore, the diffusive luminos- 
ity barely increases around the orbit in the model of log Tph 
(K)=3.9, in which the drag luminosity deposition is 1.8% of 
the total flux. 

Comparing these envelope models (black lines) with those 
without the companion effects (green and blue lines), we see 
that the envelope matter is significantly re-distributed due to 
the companion's effects. In the presence of the compan- 
ion, the Eddington luminosity, outside the orbit is Lndd = 
4^^cG{M^fJD+M^.omp)/K, which is larger than that without 
companion Lndd = 47rcGMwD / This difference causes dif- 
ferent energy flux as shown in Figure |2l 5.3 x 10^^ erg s"' 
in the model of black thick lines, whereas 4.3 x lO"'^ erg s"^ 
in the green thick lines. Therefore, these two envelopes are 
in different hydrostatic balance. A lower luminosity leads to 
a larger local super-Eddington region which causes a wider 
density-inversion region. The solution depicted by red lines 
has a photosphere just outside the companion orbit. Its pho- 
tospheric luminosity is close to those of the solutions with the 
companion, but the structure is something between the two 
solutions with/without the companion. 

Figure[3]compares structures of envelope models forMwo = 
0.55 Mq among companion masses of Mcomp = 0.0, 0.1, 0.2, 
0.3, and 0.4 Mq. For a larger companion mass, the density 
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Fig. 4. — Comparison of the static and optically thick wind solutions with 
different rotation effect for the envelopes on a 0.55 Mq WD with a 0.4 Mq 
MS companion. Static solutions are depicted by colored lines, from upper 
to lower, T] =0 (blue dashed line), 0.1 (green solid line), 0.3 (blue solid line), 
0.383 (red solid line), and 0.45 (purple dash-dotted line). Dotted black lines 
denote a static solution without the companion effects, which is almost over- 
lapped with the static solution of / = 0.383 (red solid line). All the static so- 
lutions have a photospheric temperature of \ogTph (K) = 4.3. Optically thick 
wind solutions are depicted by the black solid Hne with a black dot that indi- 
cates the critical point at Ra = 0.3 Rq ■ Three solutions of = 0.3, 0.383, and 
0.45 (the lower line in v, but the upper lines in T and p) are almost overlapped 
with each other as well as overlapped with the static solution of rj =0 (blue 
dashed line) inside the critical point. The orange dash-three-dotted lines de- 
note a test solution with different rotation law (same as orange lines in Figure 
1) which will be discussed in Section l53] 

inversion layer is smaller. In this way, the companion gravity 
has a role of redistribution of the envelope matter through dif- 
ferent hydrostatic balance of gravity and pressure-gradient; a 
large luminosity causes larger radiation pressure gradient so 
that the matter is pushed outward. Dotted and dashed lines in- 
dicate the models with enhanced drag luminosities which will 
be explained in Section l572l 

When the envelope rotates, the centrifugal force has an ef- 
fect to reduce the gravity by the factor /. Figure |4] shows the 
envelope structure for logTph (K)=4.3 with different angular 
velocity of 77=0.0, 0.1, 0.3, 0.383, and 0.45. For a larger 77, 
the centrifugal force becomes larger, which reduces the effec- 
tive gravity. When 77 = 0.383, the gravity of the companion 
(Mcomp = 0.4 Mq) is almost canceled as discussed earlier (see 
Figure[rii, and the envelope has a structure very close to that of 
the solution without the companion effects. If we further in- 
crease 77 to 0.45, the envelope structure does not change much 
anymore except for the outer region where the temperature 
and density quickly drops (see the purple dash-dotted lines). 

3.2. Structure of Wind Solutions 

Figure |4] also shows the structures of optically thick wind 
solutions. The optically thick winds are accelerated due to 
the peak of OPAL opacity which locates deep inside the pho- 



tosphere. The wind velocity quickly increases around the crit- 
ical point an d reaches the terminal velocity deep inside the 
photosphere dKato & Hachisu II 19941) . 

Optically thick wind solutions hardly change its structure 
once it is accelerated even if there are some changes of accel- 
eration source like opacity in the outside of the critical point. 
Figure |4] shows optically thick wind solutions for three dif- 
ferent 7], which correspond to different effective gravities as 
in Figure [T] As the critical point locates far inside the com- 
panion orbit and the winds are already accelerated there, the 
companion's gravity hardly changes the internal structure. 

3.3. Drag Energy Deposition during Nova Outbursts 

We already found that the drag energy deposition is small 
compared with the total luminosity. For example, it is only 
3.2% of the photospheric luminosity in the model of logTph 
(K) = 3.85 (black solid line in Figure [3]). This small contri- 
bution comes from the low density near the companion orbit 
because the drag luminosity is proportional to the density. Us- 
ing Equation (ITOl i. we estimate it as 

= TT X (3.3 X 10'°cm)^ X 6.5 x 10"'' g cm"^ x (410 km s"')^ 
= 1.6 X 10^* ergs"' 

= 0.030 Lph, (12) 

almost the same as our model value of 3.2%. 

This low density is caused by the opacity peak as shown 
in Figure |2l The density is almost constant around the com- 
panion orbit. Thus, the drag luminosity is almost independent 
of the orbital size. If we assume a less massive companion, 
then we have a smaller R-^. For a smaller companion mass, 
on the other hand, the companion's Roche radius is smaller, 
so the companion comes into closer to the WD and the enve- 
lope's density is higher at the companion's orbit. However, a 
smaller radius (accretion radius) effect overcomes the effect 
of slightly higher density. As a result, a 0.1 Mq compan- 
ion produces a hundred times smaller drag luminosity than 
an 0.4 Mq companion does. Thus, an 0.4 Mq companion 
gives almost an upper limit of the drag energy deposition for 
the 0.55 Mq WD model. During a nova outburst, the photo- 
spheric radius decreases with time and the density around the 
orbit slightly decreases with time as shown in Figure |2] in the 
case of static envelopes. Therefore, the drag luminosity de- 
creases with time. In the case of rotating static envelopes, the 
density around the orbit is as small as or smaller than that for 
the non-rotating envelope as shown in Figure IH so the drag 
luminosity is also too small. 

In the case of optically thick wind envelopes, the density 
around the orbit is also as small as or much smaller than those 
of the static envelopes as shown in Figure |4l During the out- 
burst, the density around the orbit decreases with time as will 
be shown later in Figure|2] In this way, we do not expect large 
drag luminosities during nova outbursts either in the static or 
in the wind evolution. 

Figure |5] shows the ratio of the drag luminosity to the pho- 
tospheric luminosity estimated for individual nova outbursts. 
A nova outburst starts from the uppermost point of each line, 
which corresponds to the peak luminosity, and moves down- 
ward as the envelope mass decreases with time. The drag 
luminosity quickly decreases with time because the density 
around the orbit quickly decreases in wind solutions (see Fig- 
ure |7]i. The black lines depict three cases of classical novae. 
V838 Her is one of the fastest classical novae, which oc- 
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Fig. 5. — Ratio of the drag luminosity to the photospheric luminosity dur- 
ing nova outbursts. The top of each curve corresponds to the optical peak. 
The ratio decreases with time. The bottom point corresponds to the stage that 
the companion emerges from the WD photosphere. The WD mass of each 
outburst model is indicated after the object name in solar mass units. Three 
black lines depict classical novae with optically thick winds: V838 Her (1.35 
Mq), V1974 Cyg (1.05 Mq), and GQ Mus (0.7 Mq). A red line is the com- 
posite model for V723 Cas (0.55 Mq), consisting of a static phase (top dotted 
part) and a wind phase (bottom solid pail). The dotted part of V723 Cas rep- 
resents the transition from the static to wind phase. Blue lines depict helium 
novae: V445 Pup (1.37 Mq) and a model for a 1.0 Mq WD. See text for 
more details. 

curred on a very massive WD, thus the envelope mass is very 
small. Here we adopt the envelope structure take n from the 
light c urve model of a 1.35 Mq WD calculated bv iKato et alj 
( 120091 model 2 in their Table 3) and a temporarily adopted 
0.45 Mq companion at the separation of a = 1.4 Rq. The 
lowest point of the line corresponds to the stage at which the 
companion emerges from the photosphere and the companion 
effects vanish after that. In the same way, we have estimated 
the drag luminosities for the moderat ely fast nova VI 974 Cyg 
(a 1.05 Mq WD model taken from iHachisu & Kato Il2006h 
with an 0.2 Mq companion at a = 0.85 Rq, and for one of the 
slow est classical novae GQ Mus (an 0.7 Mq WD model taken 
froml Hachisu et al. | |2008^ with an 0.1 Mq MS companion at 
a = 0.6 Rq. Here we adopt these parameters of the secondary, 
assuming they are in zero-age main sequence, which may not 
be accurate for each object. It is, however, enough for our pur- 
pose in estimating the drag luminosity, because it is not much 
affected by the choice of the secondary. 

V838 Her and GQ Mus are the fastest and slowest classical 
novae, respectively. Thus, many other classical novae may fall 
in between these two lines like V1974 Cyg. We realize that 
the drag energy deposition is quite small. Even in the very 
early phase of the outburst, the ratio reaches only a few to 
several percent. Therefore, we reconfirm that the drag energy 
deposition does not contribute to the luminosity in classical 
novae. 

Figure|5]also shows the case of a "transition nova," which is 
a model for V723 Cas (see the next section for more details). 
The binary consists of an 0.55 Mq WD and an 0.4 Mq MS 
companion. The contribution of drag luminosity is also very 
small throughout the outburst. 



The blue lines depict the case of helium nova outbursts. 
V445 Pup is only the know n helium nova for w hich we adopt 
a 1 .37 Mq WD model from lKato et"an ( 120081) . The compan- 
ion is assumed to be a l.OM© helium star and a = 2.3 Rq. For 
comparison, we also added another case of helium nova of a 
1 .0 Mq WD (model taken from Kato et al. 2008), which has 
no observational counterpart but it is a theoretical purpose. 
For the companion, we assume an 0.8 Mq helium star, and 
fl = 2 Rq. As shown in Figure |5] the drag luminosity is still 
very small. 

4. APPLICATION TO SLOW NOVAE 
4.1. Transition from Static to Wind Evolutions 

IKato & Hachisu I (120091) pointed out that two different types 
of nova evolutions, i.e., one is the evolution with optically 
thick winds and the other is without, can be realized in slow 
novae of a certain range of WD masses, ^ 0.5-0.7 Mq. 
For example, evolution of GQ Mus is explained by the se- 
quence of optically th ick wind solutions on a ^ 0.7 Mq WD 
(Hach isu et al. Il2008h while the evolution of PU Vul is de- 
scrib ed by a sequence o f hydrostatic solutions on a ^ 0.6 Mq 
WD (IKato et al ll20Tlh . 

A remarkable difference between these two evolutions ap- 
pears in the optical light curve (Kato 2011). The wind-type 
novae have a sharp peak in the optical light curve, because 
massive optically thick wind carries away a large part of the 
envelope mass and the nova light curves decay quickly. On the 
other hand, no-optically thick wind novae have a long-lasted 
flat optical peak before the magnitude slowly decays, because 
the nova evolves very slowly and it stays at an extended low- 
temperature stage for a long time, which makes a long-lasted 
fl at optical peak. 

IKato & Hachisu I (120091) presented an idea that a transition 
from the static evolution to the wind evolution could occur 
during an outburst. In such a case, the nova shows a flat op- 
tical peak with no indication of strong mass loss in the early 
phase of the outburst, fo llowed by the decay like in normal no- 
vae with strong winds. I Kato & Hach isu (2009) further sug- 
gested that such a transition accompanies some activities like 
oscillatory behavior in brightness, because the internal struc- 
tures of static/wind solutions are very different and the relax- 
ation process may cause some oscillatory features. 

Such a transition, however, may need some triggering 
mechanism. PU Vul has a flat optical peak that lasted for 8 yr 
This indicates that the static evolution is stable in a timescale 
of 10 yr. We suppose that a transition to a wind evolu- 
tion was not triggered in PU Vul because the structures of the 
static envelope are very different from those of the wind solu- 
tions. The static solution, however, changes its structure when 
a companion star disturbs the structure of a static envelope as 
shown in Section [3] This encourages us to further investigate 
our idea that the transition from the static to wind solution 
occurred in some slow novae. 

4.2. Comparison of V723 Cas, HR Del, and V5558 Sgr with 

PUVul 

Figure |6] shows the observational light curve of PU Vul, 
V723 Cas, HR Del, and V5558 Sgr PU Vul showed a flat 
optical maximum that lasted for 8 yr, neglecting the deep 
eclipse at f = 2.2 yr. The light curves of V723 Cas, HR 
Del, and V5558 Sgr are very simil ar to each other as al- 
ready pointed out by many authors (Friediung et al. "19921 
iFriediung & liiima . .2002; .Evans et al. . .2003; Munari et alH 
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Fig. 6. — Comparison of light curves among PU Vul, V723 Cas, HR Del, 
and V5558 Sgr PU Vul is shifted upward by 7 mag and HR Del and V5558 
Cyg downwai'd by 8 and 7 mag, respectively. The upper timescale is for 
PU Vul and lowe r one is for the other three novae, (a) PU Vul: data taken 
from K ato et aTl <2piT). The dips at ? = 2.2 yr and 15.7 yr are eclipses, (b) 
V723 Cas: diamonds (Chochol & Pribulla 1997), filled small circles (lAUC 
Nos. 6213, 6214, 6227, 6233, 6256, 6283, 6331, 6358, 6428), ope n circles 
(AAVSO,V-mag), d ots (A AVSO. visual), (c) HR Del: o pen cir c les jStokes I 
[TW?!: IPoEn rigw: Mia 1967; Onderlicka & VetesniTf* ^1963; 'O'Connelll 
[l%8:Grvgar 1969; Mollerus 1969; Mannerv 1970; Barnes & Evans 1970, 
V), dots (AAVSO, visual), (d) V5558 Sgr: open circles (AAVSO, V), dots 
(AAVSO, visual). The sohd lines indicate the composite light curve model 
of 0.55 Mq WD with X = 0.55, C+0 = 0.2, Z = 0.02 (red line), and 0.6 Mq 
WD with the solar composition (black line). The arrows indicate the switch- 
ing point from a static to a wind evolution: logT (K)=3.88 (red) and 3.93 
(black). 

120071) . These novae have multiple peaks, showing an oscil- 
latory variation around a certain magnitude. After that the 
magnitudes declined rather smoothly. 

Resemblances in spectral evolutions were also suggested. 
V723 Cas shows very narrow emission/absorption lines in the 
spectra at the pre-maximum flat peak. The spectrum gradually 
changes to a pure ab sorption F -type supergiant before the op- 
tical maximum (liii ma et al. i ri998). Note that PU Vul also 
shows F-type supergiant spectra with pure narrow absorp- 
tions, which are interpreted as static evolution (Kato et al. 
1201 II) . After the optical max imum, V723 Cas changes its 
spectrum to that of usual novae dliiima et al. 119981) wit h many 
emission lines, some of w hich show a P Cygni profile dliiima I 
120061: lEvans et al.1l2003h . After that it entered the nebular 
phase (645-677 days after the discovery). 
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The spectral evolution of V5558 Sgr before and a fter the 
maxirnum is also very similar to that of V723 Ca s dliiima I 
l2007bt iMuiiari et al. I l2007t iPoggiani I l2008l lIoToh . In the 
pre-maximum phase the nova shows a featureless F-type 
spectrum except narrow emission lines H a and H/3 and 
many weak absorption line s ( lijima 2007a ; Naito & Matsud a | 
120071; iTanaka et al. ll20lTI) . which does not resemble to those 
of classical novae. Fast hot winds (FWHM 1150 -1500 
km s~' ) appeared in the multi-peak episode sugge st that mass 
ejection is associated to the brightening ( Mu nari et al. 1120071 
Kiss & Sarneczky 2007; Tanaka et al. 2011). During the de- 
cline phase, V5558 Sgr had entered the nebul ar phase an d the 
spectra resemble to those of V723 Cas (.Poggianill2010l) . 

HR Del also showed spectral evolution unlike classical 
novae dHutchingsiri970l; ISanval1ll974l; iRafanelli & Rosin"o1 
119781) . The early spectra before maximum were characterized 
by relatively narrow emission lines with P Cyg ni absorption, 
and the spectral class of FO. Friedjung eTaT] ([l992) exam- 
ined spectral evolution before maximum and concluded that 
optically thick winds unlikely occurred before maximum and 
that the photosphere is almost stationary. HR Del showed 
several optical maxima which accompany discrete shell ejec- 
tions (RafanelH & Rosino 1978; Sanyal 1974). In the neb - 
ular phase, spectral evolution was alike to normal novae 
dRafanelli & Rosino l[l978.) . 

These spectral changes in V723 Cas, V5558 Sgr, and HR 
Del are consistent with our suggestion that these novae started 
its outburst from no optically thick wind evolution and, after 
the multi-peak phase, it changes into an optically thick wind 
evolution like usual novae. 

Table [T] summarizes the outburst properties of these four 
novae, i.e., the outburst year, duration of maximum phase, 
presence of multi-peak, information on dust formation, orbital 
period, estimated WD mass, and suggested type of the com- 
panion. For PU Vul, the maximum phase means the phase of 
flat maximum. For the other three novae, we define the 'max- 
imum phase' by the period between the end of quickly rising 
phase and the last prominent peak. We see that these 'maxi- 
mum phases' last 200-400 days, which is much shorter than 
the flat maximum phase of PU Vul (8 yr), but much longer 
than the orbital periods, so there is sufficient time for the com- 
panion to orbit many times around the WD before the enve- 
lope makes the transition. 

PU Vul is a long period binary of Poib = 4900 days and the 
companion is an M type red giant (RG). During the nova 
outburst, the WD envelope expands to 60 Rq), which i s 
much smaller than the companion orbit dKatoet al.l l20lT). 
Thus, the envelope structure is not affected by the compan- 
ion. On the other hand, V723 Cas (Porb = 0.693 days = 16.6 
hr: iGoran skii et al. 200Cj) and HR Del (Pprb = 0.214 days 
=5.14 hr: Friedjung.etaQ IMol; iKurster & Barwlg] [19881 
and references therein) are close binaries and the companion 
is a dwarf. The orbital period of V5558 Sgr is unknown, but 
suggested to be as short as normal classica l novae, because a 
large outburst amplitude (AV > 10 mag: iNaito & Matsuda I 
I2007i) suggests that the companion is not an RG. Therefore, 
in these three novae, the companion should have been deeply 
embedded in the nova envelope during the outburst. 

This is, again, very consistent with that PU Vul (a wide 
binary) underwent a qui et evolution w ithout strong winds 
throughout the outburst (.Kato et al. 11201 li) . No transition oc- 
curred in PU Vul. On the other hand, V723 Cas, HR Del, and 
V5558 Sgr (close binaries) outbursted as a quiet evolution at 
the beginning, then changed into a wind evolution. The multi- 
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TABLE 1 

Table 1. Slow Novae with Long-lasted Peak/Pre-maximum Phase 



Object Outburst Maximum Phase Multi-peak Dust" Porb ^WD*" Companion 

PUVul ... 1979 Syr no no 4900 days ~O.6M0 M giant 

V723Cas ... 1995 400 days yes no 16.6 hr 0.58,0.59 dwarf 

HRDel ... 1967 320 days yes no 5.14 hr 0.595. ~ 0.9 Mq dwarf 

V5558Sgr ... 2007 200 days yes no unknown 0.58-0.63 A/© dwarfs 

" PU Vul: there was a debate on the origin of the deep minimum, i.e., to be dust origin or eclipse, but later, it turned clear to be an eclipse. 
V723 Cas: no dust formed I Lvnch et al. I2QQ(J) . HR Del: no indication of dust in the spectral evolution iRafanelli & Rosino TT97I . V5558 
Sgr: no dust formed i Rudy et al 201 IJ). 

^ PU Vul: about 0.6 Mq by Kato et al. ( 201 1) from light curve fitting, V723 Cas: 0.58 ±0.07 Mq by lijima et al.^ :1998) from absolute 
B magnitude; 0.59 M q by H achisu & Kato ( 2004) from light curve fitting. HR Del: 0.5 95 Mq by Kurster & Bamig i 1988) from radial 
velocities; 0.9 bv lBruchln982D from radial velocities. V5558 Sgr: 0.58-0.63 Mq bv lPoggiaml f20fcl) from relation between absolute 
magnitude at maximum and the WD mass. 
^ see the text. 



pie peaks of these three nova light curves may be a relaxation 
process associated with the transition. The absence of mul- 
tiple peaks in PU Vul suggests no transition triggered, which 
may also support our idea. 

No indication of dust formation has been reported for all 
of the four novae as in Table [T] This is consistent with our 
suggestion that no/weak optically thick winds occur in these 
novae, because dust is usually formed in massive winds like 
OS And and DQ Her 

The estimated WD mass is about 0.6 Mq for all three no- 
vae, V723 Cas, HR Del, and V5558 Sgr Note that some 
works on slow novae mentioned that this 0.6 Mq is a mass 
close to the lower limit of WDs having s trong nova out- 
bursts as a therrnonuclear runaway event jFriediung et al. I 
[1991 iKato et al. Il200l iTanaka et al. ll20Tlh. Th is argument 
is based on the calculation bv lKovetz & PrialnikI ( 1985) with 
the old opacity. Recent calculations with the OPAL opac- 
ity, however, show that thermonuclear runaway triggers nova 
outbursts even in much less m assive WDs, e.g ., in 0.4 and 
0.6 M q WDs (,Shara et al. l [ T993i). 0.4 WD (lYaron et al. I 
l2005h . and 0.2 - 0.5 Mq WDs jShen etaL]l2009l) . Therefore, 
~ 0.6 Mq is not the lower mass limit for a WD having nova 
outbursts. Instead, we emphasize that ^ 0.6 Mq is a lower 
boundary for a nov a outburst having optica lly thick winds for 
solar composition (iKato & Hachisu I l2009h . Optically thick 
winds are weak, so it look s like a 'borderline classical nova' 
jFriediung & Iiiimall2002l) . 

4.3. Composite Light Curves 

Figure |6ta) shows theoretical light curves of PU Vul taken 
from Kato et al. (2011), which consist of only static solu- 
tions. The theoretical curve well represents the observational 
light curve until f ^ 11 yr, and then deviates after that when 
the nova enters the nebular phase. This theoretical curve rep- 
resents blackbody emission from the WD photosphere and 
does not include the contribution from nebular emission lines 
outside the photosphere as well as of the RG companion 
which dominates in the later phase. Therefore, the theoreti- 
cal magnitude is much lower than the observation. 

We made composite theoretical light curves which mimic 
the light curves of V723 Cas, HR Del, and V5558 Sgr Here, 
we assume that the early phase is approximated by a sequence 
of static solutions and the later phase by a sequence of opti- 
cally thick wind solutions. We made two trial light curves; 
one is the model of 0.55 Mq WD with the chemical com- 
position ofX = 0.55, Y = 0.23, CNO = 0.2, and Z = 0.02, 
and the other is an 0.6 Mq WD with the solar composition 
X = 0.7, Y = 0.28, and Z = 0.02. We switched the solution 



from static to wind when the photospheric temperature in- 
creases to log Tph (K)=3.89 for the 0.55 Mq WD or to log Tph 
(K)=3.87 for the 0.6 Mq WD model. The effect of the com- 
panion is included in the static phase assuming an 0.4 M© 
MS, but not included in the wind phase for simplicity because 
of their small effects. No rotation effects are included. 

Our trial light curves are depicted in Figure |6](b)-(d). Red 
lines indicate the model of 0.55 Mq WD and black lines the 
model of 0.6 Mq WD. We applied the above two trial light 
curves to the three novae, V723 Cas, HR Del, and V5558 Sgr. 
The switching epoch from a static to wind evolution is indi- 
cated by an arrow. Here, we assumed that the transition oc- 
curred ^ 200 days after the outburst. In the very beginning 
of outbursts, the transition unlikely occurs because the enve- 
lope structure may have not yet been close enough to that of 
a wind solution (see the red line solution in Figure |2]i. Once 
the envelope has extended to a giant size, it may take sev- 
eral dynamical timescale for the structure change from static 
to wind solution. The dynamical timescale is estimated to 
be fdyn ^ (2/GM)i/V3/2 = 21 days for r = 6 x 10'^ cm and 
M = 0.95 Mq. Note this timescale is roughly consistent with 
the magnitude variation of these novae which may associate 
with a relaxation process. Apart from these oscillatory be- 
haviors, our trial light curves are roughly consistent with the 
observed light curves of these three novae. 

Figure |2] shows envelope structures corresponding to each 
stage of the composite light curve model of the 0.55 M© WD 
(depicted by the red lines in Figure |6]l. The red solid lines in 
Figure |7] represent the static evolution in the flat maximum, 
the structure of which hardly changes during the flat maxi- 
mum, and the black solid lines corresponds to the wind evo- 
lution in the decay phase (after the arrow in Figure |6]l. The 
green dashed lines indicate the model without the compan- 
ion effects. We see that, in the presence of a companion, the 
transition from the static to wind evolution probably occurs, 
but may not occur without the companion effects because of 
a large difference between the envelope structures, that is, the 
green dashed line deviates largely from the rightmost black 
solid line. 

These two trial models adopt a parameter set of the WD 
mass, companion mass, chemical composition of the enve- 
lope, and the orbit, which may be close but not the exact val- 
ues for each nova. We see, however, our trial light curves 
reproduce a characteristic properties of observed light curves; 
the flat phase of static evolution followed by a smooth decline 
due to wind mass loss. If we choose another set of the param- 
eters, we can also reproduce the flat peak by static solutions 
and the subsequent decline phase by wind solutions with a dif- 
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Fig. 7. — Structure change in the evolution model for V723 Cas. The cor- 
responding light curve is shown in Figure |6] by the red lines. The red solid 
lines represent the static phase of the first 200 days (log Tpj, (K)=3.892 and 
logRph (cm)=12.65), during which the structure hai'dly changes. The black 
solid lines indicates the structure of optically thick wind phase in which the 
photospheric temperature increases with time from logTpi, (K)= 3.87, 4.16, 
4.37, 4.59, and 4.79. The con'esponding photospheric radius is logSpj, (cm )= 
12,61, 12.05, 1 1.65, 1 1.23, and 10.82, respectively. The dotted red Imes indi- 
cate a rotating static solution of r] = 0.2 and the green dashed lines the static 
solutions without the companion effects for the same photospheric tempera- 
ture (log Tpd (K)=3.892) as that of the red solid line. The black dashed line 

indicates the escape velocity (2GMwd/'')''^^. where Mwd = 0.55 Mq. 

ferent decay timescale. The exact parameter fitting is not the 
aim of the present work, but we expect that we will find a suit- 
able set of the parameters which reproduces the light curve as 
well as the other observational constraints. 

In the present paper, we studied nova outbursts with a ID 
approximation. The effects of a companion star are, however, 
highly aspherical, and the transition may occur in a different 
way in the different direction. This aspherical nature of the 
transition may relate to aspherical shapes of nova ejecta, but 
far beyond the scope of the present analysis. 

5. DISCUSSION 

5.1. Comparison with Other Works 

Common envelope evolution with frictional effects by a 
companion star has been studied by multi-dimensional hydro- 
dyna mic simulations in relation to th e formation of close bina- 
ries. iTaam & Bodenheimer I (119891) showed two-dimensional 
(2D) evolutions of a common envelope, consisting of a 5 M© 
RG and an engulfed 1 Mq MS star. Their results demon- 
strate that frictional process is strong enough to drive a rapid 
mass outflow from the equatorial plane. Eventually a large 
part of the envelope (namely > SM©) would be ejected and 
the companion would spiral in. A similar calculation but for 
a less massive 2 Mq RG and the same 1 Mq MS showed 
that frictional energy is large enoug h to ejec t the envelope, 
but the companion does not spiral in dTaam & BodenheimeTI 



19911). Three-dimensio nal (3D) hydrodynamic simulations by 
Ricker & Taarn] ( 120081) for a binary of a 1 .05 Mq RG and an 
0.6 Mq MS companion showed that the companion orbit is 
much reduced by frictional effects of the massive 0.69 Mq 
RG envelope. In these calculations, the companion efficiently 
interacts with the massive envelope. This can be understood 
from Equation (fTOl i. As the drag luminosity proportionally 
increases with the density, we expect very large drag lumi- 
n osities in mas s ive en velopes. 

iMacDonaidI (Il980h calculated nova outbursts on a 1.0 M© 
WD with a 1.39 x 10"'*Mq envelope and an 0.46 Mq MS 
companion. The author includes the drag luminosity in ID 
approximation, essentially the same as the present work, but 
the companion's gravity is not included and the opacity is 
Kramers' law. The author found that the expansion velocity 
reaches ~ 300 km s"', whereas 40 km s~' in the case without 
drag. The nova duration is not drastically shortened, i.e., the 
drag luminosity is not effective in mass ejection. 

The main difference between ours and MacDonald's (1980) 
is the opacity. With the OPAL opacity we have obtained the 
strong optically thick winds of high wind velocity (~ 1000 
km s~') for a 1.0 Mq WD and a short nova duration (s^ 1 
yr) as observed. Due to the low density at the o rbit, the drag 
luminosity is very small (iKato & Hachisu l[l994l) . From these 
reasons, we regard that MacDonald's (1980) results are not 
in consistent w i th our s . 

iLivio et alTI (119901) presented 2D calculations of a common 
envelope phase of a classical nova on a 1.0 Mq WD with a 
0.5 Mq companion and with an initially spherical envelope 
of mass 5 x 10~* Mq. Their calculation showed rapid mass 
outflows concentrated in the orbital plane. The mass loss rate 
reaches 1 x 10"'' M© yr"', corresponding to a timescale of 
5 yr for mass ejection, but still being much longer than the 
observed timescales of novae. 

As already p ointed out in our previous work 
(IKato & Hachisu I Il991bl) . iLivio et al. I (0990) assumed 
adiabatic gas, in which the drag energy is effectively con- 
sumed to push matter upward against the gravity, and results 
in the strong acceleration of matter If we adopt non-adiabatic 
gas, all the drag energy deposition can escape by diffusion 
process, as already shown in Section[3] Moreover, the authors 
adopted 60 x 40 grid points in which frictional energy and 
angular momentum are deposited into four zones near the 
companion, which may be insufficient to follow the envelope 
evolution. Therefore, we suggest that the adiabatic gas and 
low-resolution grid may be the reason why the authors get 
strong acceleration of matter concentrated in the equatorial 
plane while the drag luminosity is ve ry small. 

After the s e wor ks done by Mac Donald I 11980) and 
ILivio et al~l (Il990h . the OPAL opacity appeared which 
causes strong acceleration of nova envelopes. We ob- 
tained large mass-loss r ates and wind velocities as observed 
dKato & Hachisu Ill994l) . Once the optically thick winds oc- 
cur, the drag luminosity has little effects (see Section [3] and 
also Section 5 of Kato & Hachisu 1994). 

To summarize, we may conclude that the drag luminosity 
is effective in a massive envelope, but in a much less massive 
envelope, like realistic models of nova outbursts, the drag lu- 
minosity is entirely negligible. In other words, the common 
envelope evolution does not work in nova outbursts even if the 
companion is deeply embedded in the envelope. 

5.2. Spherical ( ID) Approximation 
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We assumed spherical (ID) approximation of the drag en- 
ergy deposition, in which the energy is deposited in a shell of 
radius Ro^b and thickness of IR.^. In 2D approximation, the 
energy is deposited in a torus of diameter of 2/?orb with cross 
section of ttRI. Thus, the energy deposition per soUd angle 
may be 



4TrRi^x2R, 
27r«o,-b X nRl 



= 2.8 



(13) 



times larger than in our ID approximation for the case of 0.55 
Mq WD and 0.4 Mq MS companion. Here we use the so- 
lution of logTph (K) =3.85 (the model denoted by the black 
solid lines in Figure [3]l, in which R^y^, = 7.4 x 10'" cm and 
7?a = 3.3 X 10'*' cm. In order to simulate this aspherical effect, 
we calculated a solution with enhanced drag luminosity by a 
factor of 2.8. That is, we multiply the right-hand side of Equa- 
tion (|9|l by the factor of 2.8. The resultant solution is shown 
in Figure|5]by dashed lines. The increase of the drag luminos- 
ity causes increase of the radiation pressure, which results in 
the density decrease around the orbit. Therefore, the resultant 
drag luminosity is 3.1 x lO-'^ erg s"', only 1.9 times the orig- 
inal one, 1.7 X 10^^ erg s"'. Even if we further increase this 
factor to 10 as an extreme case, the drag luminosity deposition 
increases only by 3.4 times of the original one, to 5.7 x lO-'^ 
erg s"'. With such enhanced energy deposition, the envelope 
structure changes, but this effects are relatively small, i.e., less 
than the change of the companion mass of 0. 1 Mq as shown in 
Figure [3] In actual multidimensional case, the energy deposi- 
tion due to frictional effects may escape by diffusion process 
into the pole (rotation-axis) directions, which makes the ef- 
fects further smaller. Therefore, our ID approximation may 
not be far from the multi-dimensional approximations. 

We also assumed that the companion mass is distributed in 
a spherical shell around the orbit. In the extended envelope, 
the photospheric radius is several tens times larger than the 
companion orbit as seen in Figure |3] Thus, the ID approxi- 
mation may not be so bad in the outer part of the envelope. In 
the vicinity of the orbit, 2D approximation, this may be close 
to a time average of a 3D case, gives a stronger gravity which 
may be roughly estimated as 



AttRI^, X 27;eomp ^ 3 g 



IttR, 



,1-b X 7r/?2^,^p 



(14) 



where 7?comp = 2.6 x 10 cm. Therefore, the companion grav- 
ity is much effective than our ID case. The companion masses 
of V723 Cas, HR Del, and V5558 Sgr are not known, but if 
it is as low as ^ 0.1 -0.2 Mq, as in many classical novae, 
the enhanced gravity may roughly correspond to our massive 
model of 0.4 Mq . Thus, envelope structure may not be largely 
different from that of our ID model. 

It should be however noted that the center of rotation is the 
center of mass, which is different from the position of the WD. 
In this sense, the center of our ID approximation (i.e., r = 0) 
switches from the WD center in the inner part to the center 
of mass in the outer part. The intermediate region between 
them is an interacting region with the companion, which is 
highly aspherical. Therefore, our ID model is a rather crude 
approximation to this intermediate region, which is not easily 
justified by our ID model. We need a full 3-D hydrodynamic 
simulation including radiative diffusion to obtain a definite 
conclusion to the transition. 



5.3. Rotation Law 

We also assumed a simple form of angular momentum dis- 
tribution in the extended envelope. Here we examine how this 
approximation affects the envelope structure, assuming a dif- 
ferent rotation law. Figure[T]shows a trial angular velocity dis- 
tribution (by the orange solid line) that represents co-rotation 
with the WD, i.e., a; = il, in the inner part of the compan- 
ion (r < /?orb -^comp, 1-6., log r (cm ) < 10.68), and specific 
angular momentum conservation in the outside of the binary 
(r > 3/?orb, i-C-, logr (cm ) > 11.35). In the intermediate re- 
gion we connect these two laws in an arbitrary form as shown 
in the figure. If the nova envelope expands very slowly, it is 
spun up by the interaction with the companion and gets angu- 
lar momentum ^ ]_2-1.7R^Q pei^unit ma ss (see num erical 
simulations, e.g. ISawada et al.1fl984l: Uahanara et al. l [2005). 
This may be the maximum value for actual nova outbursts, 
because expanding speed is not infinitely small (in V723 Cas, 
the rising timescale is one week). If no further acceleration 
occurs outside, the angular velocity decreases due to the spe- 
cific angular momentum conservation. We here assume that 
the angular velocity decreases according to the specific angu- 
lar momentum conservation, / = r^uj = 1 .lRl^^J^l at r > 3/?orb- 

Figure [T] also shows effective mass /Mi corresponding 
to the trial angular momentum distribution (the dash-three- 
dotted line). This effective mass /Mr is constant (= M,) at 
logr (cm ) < 10.2 as in the other models because the second 
term of Equation (O is small, and we have / = 1 . This means 
that the angular momentum at the WD surface is so small, 
that the results are independent of the initial rotation law of 
the WD envelope. On the other hand, the effective mass /Mr 
varies corresponding to the change of uj in the middle part 
and quickly rises as u! decreases in the outer part (log r (cm) 
> 11.3). 

Figure|4]shows the envelope structure corresponding to the 
trial rotation law by the orange dash-three-dotted lines. The 
distributions of the temperature and density in the inner part 
are very close to those of rapid rotation (rj = 0.383 and 0.45: 
red solid and purple dash-dotted lines). As explained pre- 
viously, a reduced effective gravity due to centrifugal force 
causes smaller Eddington luminosity in the outer part, and 
thus, a smaller diffusive luminosity. Therefore, the enve- 
lope mass redistributes to balance with a wider local super- 
Eddington region. As a result, the envelope structure becomes 
closer to the one without companion, although the structure is 
very different in the outer part due to irregular variation of 

/Mr. 

From these results, we say that our rotation law defined in 
Section lzTl mav be too simple but still useful in a qualitative 
study of possible transition from static to wind evolution of 
slow novae. 

6. SUMMARY 

Our main results are summarized as follows. 

1 . We present an idea of a transition from static evolution to 
optically thick wind evolution in low-mass WDs (~ 0.6 Mq). 
In close binaries, if the effects of a companion are included, 
the structure of a static envelope becomes close to that of a 
wind solution. This makes it easy to trigger the transition. In 
wide binaries, the effects of a companion are not important, 
and the structure is very different from that of the wind solu- 
tion. This makes the transition difficult. 

2. A transition from static to wind evolution could occur in 
slow novae such as V723 Cas, HR Del, and V5558 Sgr. The 
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transition can explain the spectral evolutions of these novae as 
well as their characteristic light curves. Violent activities like 
oscillatory behaviors in optical light curves are regarded as 
some relaxation processes associated with the transition. The 
presence of a companion deep inside the nova envelopes trig- 
gers this transition, which is consistent with the orbital sepa- 
rations in short-period binaries like V723 Cas. On the other 
hand, no transition occurs in long-period binaries like PU Vul 
because of no presence of a companion in the nova envelopes. 

3. We have calculated the drag luminosity for individual 
novae/helium novae with various WD masses. The drag lu- 



minosity is negligibly small during the outbursts. 
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